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Abstract 

(N ■ 

For space-times with two spacelike isometries, we present infinite hierarchies of exact 
t^J- ' solutions of the Einstein and Einstein-Maxwell equations as represented by their Ernst po- 

C*") . tentials. This hierarchy contains three arbitrary rational functions of an auxiliary complex 

parameter. They are constructed using the so called "monodromy transform" approach and 
our new method for the solution of the linear singular integral equation form of the reduced 
Einstein equations. The solutions presented, which describe inhomogeneous cosmological 
models or gravitational and electromagnetic waves and their interactions, include a number 
""^ ' of important known solutions as particular cases. 

cr 

5h . Introduction 

A number of solution-generating techniques are known which provide tools for the construction 
of vacuum and electrovacuum solutions of Einstein's equations for space-times with symmetries. 
These methods are based on the integrability of the symmetry reduced Einstein equations (viz. 
the Ernst equations). However, most of them were primarily designed to construct exact sta- 
tionary axisymmetric solutions for which an additional regularity condition should be satisfied 
on the axis. This condition does not apply to interacting waves or cosmological models as 
considered here. 

Apart from the completely linearizable subcase of EinsteimRosen vacuum gravitational 
waves, the only techniques which provide nontrivial tools for the construction of solutions for 
the dynamical case are the vacuum Belinskii-Zakharov inverse-scattering method the so 
called "monodromy transform" approach 0, ||, f|], and the group-theoretical approach recently 
developed by Hauser and Ernst ||. in particular, the methods of || enable the construction 
of soliton perturbations of homogeneous cosmological models and some specific solutions for 
wave interaction regions. For example, the Khan-Penrose or Nutku-Halil (|| solutions for 
the interaction region for colliding impulsive gravitational waves on a Minkowski background 
formally turn out to be two-soliton solutions on a symmetric Kasner background. 

Here we consider the monodromy transform approach and the linear singular integral equa- 
tions which arise in this context as an alternative form of the reduced Einstein equations. We 
present a new method for the solution of these equations which gives rise to infinite hierarchies 
of exact solutions. Among many other solutions, these include the particular cases mentioned 
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above together with other soliton solutions on the symmetric Kasner background and their 
non-soliton extensions. 



Integral equation form of reduced Einstein equations 

According to methods developed in 0, J3], ^j, any solution of the Ernst equations can be con- 
structed from the solution of the linear singular integral equation 

- / j^nr, C)<P (£, V, d( = -k(r) (1) 

considered here for the hyperbolic case only. The parameters £, r] are two real null space- 
time coordinates, e.g. = (x + t, x — t). These coordinates span some local region in 
the neighbourhood of some initial regular space-time point Pq: £ = £o> f] = t]o, in which local 
solutions of the reduced Einstein equations are considered. 

The integration in (|l|) is performed along the path L on the spectral plane w which consists 
of two disconnected parts L + and L_. In the hyperbolic case, these are chosen as the segments 
of the real axis in the w-plane, which go from w = £o to w = £, and from w = t]q to w = r] 
respectively. (We choose £o Vo and take £ and r\ sufficiently close to £o and % that the segments 
L ± do not overlap.) 

The integral in (Q) splits into two, one of which possesses a singular kernel of Cauchy type 
and should be understood as a Cauchy principal value integral. The integration parameter £ 
and a parameter r span both of the contours L + and L_. Sometimes it will be convenient to 
introduce suffices: C+i T + £ L + and C-i r - £ L_. 

In the integrand in (jl]), [A]^ = |(Ai e ft — A r i g ht)- This represents the jump on the contour, 
i.e. half of the difference between left and right limit values at the point £ £ L + or £ G L_ of 
some "standard" function X(^,rj,w). This function is a product of two functions A(£, r),w) = 
\+(£,w)\_(ri,w) given by 



A + = . /Jizi, A_ = J^L, (2) 
V w - 4o V w 

with the additional conditions X + \ w=OQ = X_\ w=OQ = 1. Each of these functions is an analytic 
function on the whole spectral plane w apart from the cut L + or L_ respectively, whose endpoints 
are the branching points of the corresponding function. 

In the equations (|l|), the three-dimensional complex vector function cp (£, 77, £) is unknown, 
and the right hand side k(r) is a three-dimensional complex vector function of the spectral 
parameter which may be taken to be 

kH = {i,uH,vH}, (3) 

where u(w) and v(w) are arbitrary functions. The kernel of the integral in (|j) is a scalar function 
H(r, C) given by 

H(r, = 1 + i(C - A))(u(t) - u+(C)) + agu(r)ut(C) 

-4(C-&)(C-%)v(r)vt(C) (4) 



where the dagger denotes complex conjugation: e.g. w(w) = u(w). The additional constants 
in (|) are a = (Co - m)/ 2 and (3 = (Co + %)/2. 

It is important to emphasize that the integral equations (|]), and hence the functions u(w), 
v(w) and <p (C, r/, w), only need to be evaluated on the two cuts L + and L_ in the spectral plane. 
Thus all the above vector and scalar functions of the spectral parameter are actually determined 
by pairs of functions which represent their values on these contours. For convenience we shall 
denote the values of these functions on L± by the corresponding suffices: 

fu(w)v(w)\ - { { U +(H V +H}> weL + 

iU[w) > V[wn -\{u_(w),v_(w)}, weL. [b) 



Thus, in (|T|) written in a more explicit form, we actually have two unknown vector functions 
cp ± . For any of these suffixed functions we can use also an alternative definition, for example, 

<P (£)f7>7±) = cp ± (£,?7,r). 

Using this notation, it is convenient to split the integral in ([!]) into separate integrals over L + 
and L_ and to consider separately the cases r = r + G L + and r = r_ 6 L_. It is also convenient 
to denote the four scalar kernels which appear in the integrands of (|lj) in the form 

H(t + , c + ) = H ++ (t, C), H(t + , c_) = C), 

H(t_,( + ) = W-+(r,C), W(r_,C-) = W-(r,C) 

where the functions H ++ (t, (), TC + _(t, (), TC_ + (t, () and W (r, £) can be determined explicitly 

in terms of the four functions u±(w) and v±(w) using (Q). 

To conclude our description of the structure of the master integral equations, we recall 
that the four functions u±(w) and v±(w) appearing in (|3|) and (|5|) play a significant role in 
the entire construction. They determine completely the coefficients of the integral equations 
in the electrovacuum case. In the vacuum case there are only two such functions u ± (w), as 
v±(w) = 0. As shown in ||, they characterize unambiguously every individual solution of the 
Ernst equations. Moreover, the singular integral equations ([j]) possess a unique solution for any 
given choice of analytical functions u ± (w) and v ± (w). 

We recall now also, that the general local solution of the hyperbolic Ernst equations can be 
expressed by quadratures in terms of the solution of (|l]) 

S = -l- I f [A] f [l-i(C - /3o)uHO]<P [u] &V,Od( 

IT Jl 

* = - / [\] c [l-KC- (3 W(OW [v] (t;,v,()dC (6) 

where cp ^ and <p ^, in some association with the definition (^), denote respectively the second 
and third components of the vector solutions cp of the master integral equation (|]), corresponding 
to a given choice of the monodromy data functions u±(w) and v ± (w). In a more explicit form, 
each of the the integrals in (||) should be split into two integrals evaluated over L + and L_. 



New hierarchies of solutions 

Here we will construct a class of hyperbolic solutions that is determined by the rational mon- 
odromy data 

/ x U ± (w) V ± {w) 
Q±[w) Q±{w) 

where U + (w), V + (w), Q + (w) and U_(w), V_(w), Q^{w) are arbitrary complex polynomials, 
provided u + (u>), v+(w) and u_(u>), v_(iu) do not have poles on L + and L_ respectively. 

For what follows, it is convenient to calculate some auxiliary polynomials of two variables - 
we introduce the four polynomials P±±(t, Q defined by the relations 

74±(t,C)= P±dT \° , (8) 

and four polynomials R±±(t, Q defined from them by 

^(r.o- ^ (Ti0 { :^ (c '° . (») 



In these definitions there are two sets of suffices, denoted as dotted and undotted, which should 
each be taken to be the same. Finally, it is convenient to introduce a redefinition of the unknown 
functions 



(P+(C) = - P ++ (C,C) V+(C) ' 

<p-(C) = — p <p-(0- (10) 

Hereafter we do not show explicitly the arguments £ and r\ of <p ± and A or the suffices ± at the 
points C and T ; unless it is necessary. 

A direct substitution of ([?]) into equations (||) with the use of (ffH(|i"0|) leads to the following 
convenient form of linear equations with polynomial right hand sides 



Co So 



»7 

7T2 J 



'/() 



■R+_(r + ,C) ~ ,> 




(11) 



m J Q - t_ m J -P--(C) 

»7o % 



1 /[A + ]c VH ) v + (0<*C + f ?£, 

w p ++ (c,o \ v _ iT) 



if we impose constraints on the coefficients of the rational functions (0) such that 

p + _(C,C) = P- + (C,C) = o. (12) 



This leads to a large class of explicit solutions <p ±(£,i],t) of (11) that are regular on the cuts 
L±. However, the solution of the Ernst equations needs the solutions cp ± (£, rj, t) of (|l|) to be 
regular on the cuts L±. Fortunately, all additional singularities (poles) of cp + 77, r) on L + and 
<p _(£, 77, r) on L_, which arise from the denominators in (|To|), can be avoided by the additional 
restrictions that u + (r/o) = —i/o-o and u_(£ ) = i/«o- We therefore specify 

. i C + (u>) 

u+W = -^ + ( "- r?o) ^H 

!!_(„,) = -L + ( u ,_ & )glM (13) 

a Q-[w) 

where C+(w), C-(w), Q + (w) and Q_(u>) are arbitrary polynomials. With these, the ansatz 
([l2] ) leads to the constraint C-(w) = B(w)C+(w) — 4iA(w)V+(w) and, for the polynomials in 
@, the general solution of ( |i"2|) reads 

Cf+H = -— Q+(w) + {w- Vo)C+(w) 

U_(w) = B(w) (j^QlH + ( w ~ A>)C+H) - 4i(w - t )A(w)Vl(w) 
V-(w) = A(w) (q\.(w) - ialC ] + {w)) 

Q-H = B(w) (q\(w) - iagC^M) (14) 



where the polynomials A(w), B(w), C + (w), V + (w) and Q + (w) can be chosen arbitrarily, pro- 
vided the corresponding functions u±(w), v±(w) have no poles on the cuts L + and L_ respec- 
tively. The vacuum case, which occurs when A(w) = V+(w) = and B(w) = 1, yields simpler 
expressions which involve just two arbitrary polynomials C+(w) and Q + (w). 

Returning to (|TT|), we note that the integral operators in the left hand sides can be inverted 
using the Poincare-Bertrand formula || for singular integrals 

- [ -^-cp (0 dC = f(r) & cp (r) = - [ ^C/(C) dC (15) 

This can be applied to the integrals over L + (using A + ), or over L_ (using A_). 

Since the right hand sides of ([ll]) are polynomials in r, the inversion (|l5| ) leads to the solution 
in the form 

N ± I Qk± \ 

<P±( r ) = ^2\Uk± \ Z k± (r) (16) 
fc =o V v k± J 

where N + and N_ are the maxima of the degrees of the polynomials U + , V + , Q + and U_, VI, 
Q_ respectively, Uk±, Vk±, qk± are unknown r-independent functions of £ and 77, and Zk±{r) are 
"standard" functions given by 

1 f [AJ 1 ]^ 



Z k± (r) = - ^^CdC (17) 
m Jl ± Q — r 

All these functions (integrals) can be evaluated as the residues of their integrands at £ = 00 are 
polynomials in r of degree k. 

We note now, that the vector integral equations (0) decouple into three pairs of equations 
- one pair for each of the three components of <p + and the corresponding component of cp_. 
All these pairs of equations possess the same kernels but different right hand sides. Therefore, 
substituting the expressions dH) into (0) and using (^) with ([l2|), we get three decoupled 
algebraic systems, each of order (N + 2) x (iV + 2) where = N + + N_ and for the sets of 
unknowns q k ±, u k ±, v k± respectively. However, in view of @, we need the solutions of two of 
these systems only: 

N+l 

^2 T^abUb = u A , 



T, V AB v B = v A 

B=0 

where the indices A, B = 0, 1, . . . N + 1. The column vectors u A , v A (shown below as rows) are 
composed of the coefficients of the polynomials ?7±(C) and V ± (£): 

u A = {u 0+ ,u u , . . . ,u N+ ,u _,ui_, . . . ,u N _} 

v A = {v 0+ ,v 1+ , . . . ,v N+ ,v _,v 1 _, . . .,v N _}. (19) 

Similarly, we combine the coefficients u k ±, v k± in ( |l6| ) to form the column vectors (rows) 

Ua(€,t]) = {uo+,ui + , . . . ,U N+ ,U _,Ul_, . . .,u N _} 

v A (€,V) = {vq+,v 1+ , . . . ,v N+ ,v _,vi_, . . . ,v N _} (20) 

The matrix [|X>|| consists of the blocks D ++ , D + _, D_ + , D of orders (iV + + 1) x {N + + 1), 

(N + + 1) x (N_ + 1), (N_ + 1) x (N + + 1) and (N_ + 1) x (N_ + 1) respectively. Their components 
are determined by the integrals: 



TTl J P++{C,Q 



£0 



VO 

(D-MO = - [ [Kk ^l^} Zl + (0 d( (21) 
m J P++{(,() 



5o 



(D__) kl (rj) = 5 kl + -.l [\-}c { u~~)^ Z l-(OdC 



TTiJ 1 P__(C,C) 

where (R ±± )k are the coefficients in the expansions R + ±(t, () = J2k=o(R+±)k(() Tk an d R-±( T > = 

To calculate the final expressions for the Ernst potentials, we need to evaluate the additional 
sets of integrals 

Jk+{0 ~ Vi J [A+k R^O fc+(0 C 

Co 

juv) - l - /[A-ic q - (c) -; K -$ )vlK) zuo*, 

Vo 

and to combine them into one row vector 

J A = {Jo+, J\+, • • • , Jn+i Jo-i J1-1 • • • 1 Jn-}- (22) 
Let us also define two additional (N + 2) x (N + 2) matrices 

Gab = Dab ~ ^iu^Jg, 3~ ab = Dab ~ ^IvaJb- (23) 

All integrals determining the components of the matrices Gab, Fab and T> AB can be evaluated in 
terms of the residues of their integrands at the zeros of P ++ (w, uu) and P__ (w, w) and at w = 00. 
We then have 

_ det \\Gab\\ $ = det \\T AB \\ ^ 

det||X> AB ||' det||P AS ||' 1 ' 

which are the final expressions for the Ernst potentials. These solutions generally possess es- 
sentially nonlinear properties. They are not trivial time-dependent analogues of any stationary 
axisymmetric solutions with regular axis of symmetry which have different structures of mon- 
odromy data. The expressions (24) generally are not rational functions of £, rj. 

When evaluating explicit examples, it may be noted that solutions with a diagonal metric 
occur when u± = — u ± . The plane symmetric (type D) Kasner metric with £ = —a/aQ is 
obtained using the constants u + = —i/a>o, u_ = i/ao and v ± = 0. The Khan-Penrose solution 
@ for colliding plane impulsive gravitational waves is obtained with v+(w) = v_(u;) = and 

w-a + w-a- . . 

u.+ {w) = ik + — , u_(w;) = iK- — (25) 

w — b + w — b- 

when the constants a±, b± and k± are real. The nondiagonal Nutku-Halil solution S for non- 
colinear impulsive waves is obtained from the same expression using complex constants. This 
explicitly demonstrates that the above method is applicable to both the linear and nonlinear 
cases. 
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